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The  addition  of  a  fast  auxiliary  power  source  like  a  supercapacitor  bank  in  fuel  cell-based  vehicles  has 
a  great  potential  because  permits  a  significant  reduction  of  the  hydrogen  consumption  and  an  improve¬ 
ment  of  the  vehicle  efficiency.  The  Energy  Management  Strategies,  commanding  the  power  split  between 
the  power  sources  in  the  hybrid  arrangement  to  fulfil  the  power  requirement,  perform  a  fundamental  role 
to  achieve  this  objective.  In  this  work,  three  strategies  based  on  the  knowledge  of  the  fuel  cell  efficiency 
map  are  proposed.  These  strategies  are  attractive  due  to  the  relative  simplicity  of  the  real  time  imple¬ 
mentation  and  the  good  performance.  The  strategies  are  tested  both  in  a  simulation  environment  and  in 
an  experimental  setup  using  a  1.2-kW  PEM  fuel  cell.  The  results,  in  terms  of  hydrogen  consumption,  are 
compared  with  an  optimal  case,  which  is  assessed  trough  an  advantageous  technique  also  introduced 
in  this  work  and  with  a  pure  fuel  cell  vehicle  as  well.  This  comparative  reveals  high  efficiency  and  good 
performance,  allowing  to  save  up  to  26%  of  hydrogen  in  urban  scenarios. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Given  the  serious  environmental  problems  and  the  anticipated 
fuel  shortage  for  the  next  decades  it  is  important  to  find  more 
efficient  forms  of  using  the  power  resources,  affecting  minimally 
the  environment.  According  to  different  sources,  such  as  [1  ],  global 
oil  reserves  are  only  sufficient  for  around  40  years  with  the  cur¬ 
rent  level  of  oil  production.  However,  the  global  economy  is  rapidly 
increasing  with  a  subsequent  increase  in  oil  consumption:  the  total 
primary  energy  consumption  has  grown  from  6128  Mtoe  in  1973  to 
11435  Mtoe  in  2005;  meanwhile,  the  C02  emission  has  grown  from 
15661  Mt  to  7136  Mt  in  the  same  period  [2].  Thus,  oil  reserves  will 
therefore  be  probably  exhausted  in  a  much  shorter  time  period. 

The  main  advantage  of  fuel  cell-based  vehicles  as  opposed  to 
those  with  internal  combustion  engine  is  that  in  the  former  the 
energy  conversion  is  direct  (without  combustion),  allowing  a  higher 
efficiency  and  practically  null  emission  of  polluting  agents.  The 
efficiency  of  a  fuel  cell-based  vehicle  using  hydrogen  is  twice  the 
efficiency  of  a  gasoline  vehicle  [3-4].  Moreover,  unlike  the  internal 
combustion  engines,  the  fuel  cell  efficiency  is  also  high  at  partial 
loads:  a  vehicle  running  in  urban  or  suburban  areas  is  demanding  a 
small  fraction  of  the  rated  fuel  cell  power  most  of  the  time  [5].  Thus, 
fuel  cell-based  vehicles  have  long  term  potential  to  be  the  main¬ 
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stream  vehicle  in  the  future  because  they  have  almost  zero  emission 
and  compatible  driving  range  to  internal  combustion  vehicles  [6]. 

Fuel  Cell  Hybrid  Vehicles  (FCHVs)  is  a  promising  fuel  cell  appli¬ 
cation  that  has  taken  more  and  more  importance  in  the  last  years 
and  is  considered  the  most  attractive  long-term  option  for  pas¬ 
senger  cars  [7].  Hybridization  in  fuel  cell-based  vehicles  consists 
in  adding  a  supplementary  energy  storage  element  (e.g.,  a  bat¬ 
tery  or  a  supercapacitor  bank)  to  the  primary  power  source  (i.e., 
the  fuel  cell  system)  in  order  to  adequate  in  an  optimum  way 
the  energy  generation  to  the  consumption.  This  procedure  has 
important  advantages,  allowing  a  greater  reduction  of  the  hydrogen 
consumption.  Depending  on  how  the  necessary  power  is  obtained, 
a  minimization  of  the  hydrogen  consumption  can  be  achieved:  the 
Energy  Management  Strategies  (EMSs)  are  algorithms  which  deter¬ 
mine  at  each  sampling  time  the  power  generation  split  between  the 
fuel  cell  system  (FCS)  and  the  energy  storage  system  (ESS)  in  order 
to  fulfil  the  power  balance  between  the  load  power  and  the  power 
source. 

The  EMSs  in  a  FCHV  is  an  open  issue,  existing  different 
approaches  in  the  literature.  In  the  work  of  Thounthong  et  al.  [8],  the 
strategy  is  based  on  DC  link  voltage  regulation  though  the  control  of 
the  power  converters  but  does  not  deal  with  the  system  efficiency. 
Jung  et  al.  [9]  propose  a  control  strategy  with  two  objectives,  to 
obtain  high  efficiency  and  to  maintain  a  minimum  state  of  charge  in 
the  batteries,  giving  priority  to  the  second.  On  the  other  hand,  other 
literature  works  are  focused  on  strategies  based  on  optimization 
theory.  Paganelli  et  al.  [10]  present  a  method  based  on  a  con- 
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Nomenclature 

APfcsjaii  rate  maximum  fall  rate  of  PfCS 

Apprise  rate  maximum  rise  rate  of  PfCS 

AT 

sampling  time 

m 

efficiency  of  the  Boost  converter 

Pb/b 

efficiency  of  the  Buck/Boost  converter 

PI  ess 

efficiency  of  the  ESS 

Pfcs.max 

maximum  efficiency  of  the  FCS 

ConsH2 

hydrogen  consumption 

Cr 

capacitance  of  the  supercapacitors 

DC 

direct  current 

Pcap 

maximum  energy  that  the  ESS  is  capable  to  store 

EMS 

Energy  Management  Strategy 

ESS 

energy  storage  system 

FC 

fuel  cell 

FCS 

fuel  cell  system 

FCHV 

Fuel  Cell  Hybrid  Vehicle 

kess 

constant  of  the  ESS 

lhvH2 

low  heating  value  of  hydrogen 

M 

mass 

Nc 

number  of  sampling  time  in  a  driving  cycle 

P 

power 

PfCS, hi 

parameter  in  the  EMSs 

Pfcs.low 

parameter  in  the  EMSs 

Pfcs.max 

rated  fuel  cell  power 

Rd 

internal  resistance  of  the  supercapacitors 

SoE 

state  of  energy  in  the  ESS 

tc 

duration  of  a  driving  cycle 

trol  strategy  called  Equivalent  Consumption  Minimization  Strategy. 
This  strategy  demonstrates  to  be  robust  under  a  wide  range  of 
operating  conditions.  Nevertheless,  the  efficiency  is  sensible  to 
several  parameters.  In  the  same  line,  Rodatz  et  al.  [11]  imple¬ 
ment  a  real  time  control  by  using  the  same  concept  of  equivalent 
hydrogen  consumption,  showing  experimental  results.  Another 
optimization  algorithm,  presented  in  [12],  diminishes  a  cost  func¬ 
tion  depending  only  on  system  variables  considered  at  the  present 
time  but  the  conclusion  is  that  the  efficiency  is  not  significantly 
improved. 

With  regard  to  the  energy  storage  element  in  a  FCHV,  some  pub¬ 
lications  in  the  literature  analyze  the  feasibility  and  performance  of 
supercapacitors.  Karden  et  al.  [  13  ],  analyzing  energy  storage  devices 
for  hybrid  electric  vehicles,  conclude  that  supercapacitors  are  effec¬ 
tive  for  pulse  power  applications.  Kotz  et  al.  [14]  analyze  a  pilot 
system  on  small  scale  that  represents  a  hybrid  vehicle  with  fuel  cell 
(FC)  and  supercapacitors,  demonstrating  that  supercapacitors  are 
excellent  energy  storage  devices,  which  can  be  efficiently  used  in 
vehicle  applications  for  energy  braking  recovering  and  peak  power 
reinforcing.  Another  approach  using  supercapacitors  is  found  in  the 
work  of  Caux  et  al.  [15,16]  in  a  tramway  application.  Schiffer  et  al. 
[17]  and  Kim  and  Peng  [18]  also  show  the  benefits  of  using  a  fuel 
cell-supercapacitor  system  for  FCHVs. 

The  main  objective  of  EMSs  is  to  improve  the  efficiency  reduc¬ 
ing  the  hydrogen  consumption  along  the  driving  cycle.  To  find  a 
global  optimal  solution,  control  techniques  where  a  minimization 
problem  is  resolved  have  been  studied,  e.g.  in  [19].  In  general,  these 
techniques  do  not  offer  a  causal  solution,  and,  in  consequence,  are 
not  feasible  because  it  is  assumed  that  the  future  driving  cycle  is 
entirely  known  a  priori.  Nevertheless,  their  results  can  be  used  to 
make  comparative  analysis  of  the  performance  of  new  strategies  in 
study.  On  the  other  hand,  strategies  which  deal  with  local  optimiza¬ 
tion  are  convenient  for  real  implementation.  Another  approach, 
particularly  convenient  to  work  in  real  time  operation,  is  the  type 


of  strategies  where  the  system  operation  is  based  on  heuristic  rules 
derived  from  the  knowledge  of  the  performance  of  the  system’s 
components.  These  strategies  can  achieve  a  comparable  perfor¬ 
mance  to  those  strategies  where  some  optimization  technique  is 
utilized. 

In  this  work,  three  Energy  Management  Strategies  for  FCHVs 
are  presented.  The  strategies  are  based  on  the  knowledge  of  the 
fuel  cell  efficiency  map:  two  of  them  are  heuristic  type  strategies 
and  the  third  strategy  is  based  on  constrained  nonlinear  program¬ 
ming.  Besides,  a  novel  method  to  asses  the  optimal  consumption 
knowing  the  future  driving  profile  is  presented.  This  last  result  is 
used  to  compare  the  performance  of  the  proposed  strategies.  The 
FCHV  in  study  is  based  on  a  small  car  with  the  energy  storage  sys¬ 
tem  composed  by  a  supercapacitor  bank  and  the  entire  system  is 
modelled  using  ADVISOR  (ADvanced  Vehicle  SimulatOR  is  a  tool¬ 
box  developed  by  the  National  Renewable  Energy  Laboratory  with 
the  aim  of  analyzing  the  performance  and  fuel  economy  of  conven¬ 
tional,  electric,  and  hybrid  vehicles  [20,21  ])  in  MATLAB/SIMULINK. 
The  proposed  strategies  are  tested  both  in  simulation  and  in  an 
experimental  setup,  showing  an  auspicious  performance. 

The  organization  of  the  paper  is  as  follows.  In  Section  2,  Fuel 
Cell  Hybrid  Vehicles  are  introduced  showing  the  advantages  and 
the  mechanism  that  allows  improving  the  efficiency.  Then,  the  pro¬ 
posed  Energy  Management  Strategies  are  introduced  in  Section  3. 
In  Section  4,  the  simulation  results  obtained  using  these  strate¬ 
gies  are  presented,  including  the  description  of  the  procedure  to 
asses  the  minimum  consumption.  In  Section  5,  the  experimen¬ 
tal  validation  is  presented,  describing  the  experimental  setup  and 
showing  the  comparative  results.  Finally,  the  discussion  of  the 
results  and  the  conclusions  of  the  work  are  in  Sections  6  and  7, 
respectively. 

2.  Fuel  Cell  Hybrid  Vehicles 

In  a  Fuel  Cell  Hybrid  Vehicle  (FCHV),  the  powertrain  is  composed 
by  the  fuel  cell  system  and  an  ESS  in  order  to  increase  the  efficiency 
of  the  integrated  system  (Fig.  1 ).  Given  a  certain  load  power  PioadW * 
this  can  be  supplied  with  some  power  from  the  fuel  cell  system, 
Pfcs(t),  being  the  rest  of  power  supplied  by  the  energy  storage  sys¬ 
tem,  Pess(t).  The  power  balance  in  the  DC  bus  must  be  fulfilled  at 
every  time 

Pfcs(t)m  +  Pess(t)riB/Br)ess  =  P/oadiO  Vt,  (1 ) 

where  rjB,  rjB/B,  and  rjess  are  the  efficiency  of  the  power  converter  that 
connects  the  fuel  cell  with  the  bus,  the  efficiency  of  the  power  con¬ 
verter  that  connects  the  ESS  with  the  bus,  and  the  efficiency  of  the 
ESS,  respectively.  The  control  of  the  DC/DC  converters  is  beyond  the 
scope  of  this  paper.  We  assume  that  the  converters  are  controlled 
and  their  efficiencies  are  known. 

In  a  FCHV,  three  different  mechanisms  allow  the  efficiency 
improvement 

1 )  In  automotive  applications,  it  is  possible  to  recover  energy  from 
regenerative  braking,  allowing  a  considerable  improvement  in 
the  hydrogen  economy.  The  recovered  energy  from  regenerative 
braking  may  represent  an  important  fraction  of  the  economy 
improvement.  However,  the  amount  of  energy  that  is  possible 
to  recover  is  strongly  dependent  on  the  driving  profile.  In  urban 
driving  cycles,  in  which  there  is  considerable  braking,  the  per¬ 
centage  of  regenerated  energy  is  much  greater  than  in  non-urban 
cycles  where  this  percentage  is  quite  small. 

2)  The  efficiency  of  a  FCS  is  high  at  partial  loads.  Therefore,  it  is 
convenient  to  perform  an  energy  management  in  which  the  fuel 
cell  is  preferably  operated  in  its  maximum  efficiency  zone.  Thus, 
taking  advantage  of  an  auxiliary  power  source,  e.g.  the  storage 
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Fig.  1.  FCHV  scheme. 


system,  and  an  adequate  energy  strategy,  it  is  possible  to  main¬ 
tain  the  fuel  cell  operating  at  its  maximum  efficiency  region  most 
of  the  time  and  reduce  the  hydrogen  consumption. 

3)  In  a  pure  fuel  cell  vehicle,  the  fuel  cell  must  be  sized  to  meet  the 
largest  load.  One  important  advantage  of  hybridization  is  that 
it  allows  a  diminution  of  the  rated  power  of  the  fuel  cell  stack 
in  these  types  of  applications.  Therefore,  it  is  possible  to  meet 
the  demand  with  a  smaller  FCS  in  combination  with  an  energy 
storage  system.  Besides,  downsizing  the  FCS  has  an  impact  on 
the  hydrogen  economy  in  three  ways:  (i)  downsizing  leads  to 
improvements  in  the  hydrogen  economy  when  the  peak  effi¬ 
ciency  of  the  FC  stack  is  shifted  to  better  match  with  the  power 
requirements  for  a  given  load  cycle,  (ii)  downsizing  leads  to  a 
smaller  compressor  reducing  the  parasitic  losses,  and  (iii)  down¬ 
sizing  reduces  the  total  mass  of  the  vehicle. 

Although  the  main  motivation  for  introducing  hybridization  in 
fuel  cell  systems  is  to  improve  the  hydrogen  economy,  the  benefits 
are  not  limited  to  that.  In  fact,  hybridization  is  also  useful  to  solve 
two  important  problems  in  fuel  cell  control 

1.  The  dynamics  of  FCS  are  relatively  slow,  mainly  because  of  the 
dynamics  of  the  air  compressor  and  the  manifold-filling  dynam¬ 
ics  [22].  In  this  context,  an  energy  flow  from  an  energy  storage 
device  with  high  specific  power  output  is  useful  to  improve  the 
transient  response.  In  the  literature,  this  issue  is  not  completely 
solved.  In  [23],  it  is  concluded  that  it  is  feasible  to  use  a  FCS  in 
a  load-following  mode  without  any  storage  device.  In  that  work, 
the  conclusions  are  based  on  results  obtained  from  experimental 
data  working  with  a  scaled-down  FCS  through  the  Federal  Urban 
Driving  Schedule  (FUDS).  It  is  remarked  that  it  is  possible  to  pull 
short  bursts  of  power  from  the  stack,  due  to  the  charge  in  the 
double  layer  (CDL)  capacitance.  Flowever,  several  publications, 
for  example  [14],  [24]  or  [4],  conclude  that  the  use  of  energy 
storage  devices  is  beneficial  from  both  the  system  dynamics  and 
the  hydrogen  economy. 

2.  Flybridization  can  help  avoiding  oxygen  starvation  by  means  of 
the  energy  supply  from  the  energy  storage  elements.  Thus,  the 
peaks  of  actual  current  demand  to  the  fuel  cell  are  lower  than 
in  the  pure  fuel  cell  case  where  no  auxiliary  power  source  is 
present.  Oxygen  Starvation  is  a  complicated  phenomenon  that 
occurs  when  the  oxygen  partial  pressure  falls  below  a  critical 


level  at  any  point  in  the  cathode.  It  can  be  produced  by  an  abrupt 
increase  in  the  current  demand  (that  produces  a  rapid  increase 
in  the  oxygen  consumption  due  to  the  cathode  electrochemical 
reaction)  and  an  insufficient  air  supply.  This  phenomenon  can 
be  observed  by  a  rapid  decrease  in  the  cell  voltage  and  in  severe 
cases  can  produce  a  short  circuit  and  a  hot  spot  on  the  surface  of 
the  membrane  cell  [25].  The  phenomenon  is  not  well  understood 
and  more  information  from  electrochemical  studies  is  needed 
to  be  able  to  quantify  the  critical  level  of  oxygen  concentration 
below  which  permanent  damage  is  produced  or,  perhaps  more 
probably,  the  combination  of  current  level  variation  and  its  time 
duration. 

3.  Energy  Management  Strategies  for  FCHVs 

In  this  work,  a  search  for  new  EMSs  in  FCFIVs  is  addressed  and, 
despite  the  Fuel  Cell  Flybrid  Systems  have  other  applications,  we 
concentrate  our  attention  in  automotive  applications.  The  objective 
is  to  improve  the  performance  of  a  FC  hybrid  vehicle  in  terms  of 
hydrogen  economy  and  transient  response. 

In  Fig.  1  is  represented  the  proposed  topology  of  the  system, 
where  the  FCS  is  connected  to  the  DC  bus  through  a  step-up  power 
converter  (Boost  converter),  whereas  the  ESS  is  connected  to  the 
DC  through  a  bidirectional  power  converter  (Buck-Boost  converter) 
allowing  the  energy  flow  from  regenerative  braking.  With  regard  to 
the  load,  which  is  an  AC  induction  motor,  it  is  fed  through  a  DC-AC 
inverter.  This  topology  is  common  in  the  three  strategies  proposed. 
In  the  same  Fig.  1 ,  the  energy  flows  in  the  FCFIS  are  also  represented: 
the  regenerated  energy  from  the  load  to  the  ESS  bank,  the  charging 
energy  from  the  FCS  to  the  ESS,  the  boosting  energy  from  the  ESS 
to  the  load,  and  the  fuel  cell  energy  that  directly  supplies  the  load 
from  FCS. 

3.1.  Objectives  of  EMSs  in  FCHVs 

The  hydrogen  consumption  map  of  a  FCS  clearly  reveals  a  zone, 
below  a  limit  power,  where  the  hydrogen  consumption  is  very  high 
meanwhile  the  zone  above  this  limit  power  has  a  significant  lower 
consumption  as  can  be  seen  in  Fig.  2.  The  principal  objective  of 
the  Energy  Management  Strategy  is  to  minimize  a  cost  function ] ,  a 
mathematical  expression  that  represents  the  cumulative  hydrogen 
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consumption  during  the  time: 


min/(X)  subject  to  H(X)  =  0  and  G(X)  <  0, 
where  the  vector  X  is 


X  = 


PfcsW 
Pess{k ) 


(2) 

(3) 


This  means  that  the  EMS  has  to  determine  the  optimal  value  of 
the  fuel  cell  system  power,  PfCS(t),  and  the  energy  storage  system 
power,  Pess(l<),  in  order  to  minimize  the  cost  function  J(X) 


Nc 

J(X)  =  5>(X(k))  AT 

k= 0 


(4) 


where  Nc  =  tc/AT  is  the  duration  of  the  driving  cycle,  assuming  a 
constant  sampling  period  AT  =  1  s,  and  F(X(k))  is  the  hydrogen  con¬ 
sumption  according  to  the  hydrogen  consumption  map  of  the  FCS 


F(X(lc))  =  ConsH2(Pfcs(k))  (5) 

The  constraint  H(X)  =  0  in  (2)  means  that  at  every  time  the  power 
load  in  the  DC  bus  must  be  satisfied  by  the  power  from  the  FC  and/or 
the  ESS 


PfcsW  VB  +  Pess{k )  PIb/B  Vess  —  Pload(k)  k  —  0,  1 ,  .  .  .  ,  Nc  (6) 

where  the  efficiencies  r]B,  rjB/B,  and  rjess  are  assumed  to  be  constants. 

On  the  other  hand,  the  constraint  G(X)  <0  in  (2)  states  the  con¬ 
straints  in  the  FCS  and  the  ESS  power  for  all  k 


^  dPfcs{k ) 

APfcsJall  rate  —  ^  —  APfcs^  riserate 

Pess,dischrg  max(^°P(^))  <  Pess(k )  —  Pess,chrg  maxi 1 

SoPmin  —  SoE(l<)  <  SoEmax 


(7) 

(3) 

(9) 

(10) 


The  inequality  (7)  states  the  restrictions  in  the  FCS  power.  The 
maximum  power  is  limited  by  the  FCS  rated  power,  whereas  the 
minimum  FCS  power  command  must  be  limited  to  a  value  below 
which  it  is  not  suitable  to  operate  because  the  parasitic  load  is  too 
large,  reducing  the  system  net  power. 

There  is  a  phenomenon  of  delay  between  the  onset  of  the  load  on 
the  FCS  and  the  response  of  the  reactant  supply  system  which  may 
lead  in  a  “starvation”  of  reactants  in  the  FCS.  This  may  be  avoided 


by  restricting  the  dynamics  required  by  the  load  [12].  The  same 
approach  is  followed  in  [26]  and  [14]  where  the  maximum  posi¬ 
tive  rise  is  limited  to  avoid  damage  in  the  stack.  Therefore,  in  our 
work,  the  FCS  power  load  is  increased  no  faster  than  a  certain  power 
rise  rate  (APfcsrise rate).  Besides,  we  propose  to  operate  the  FCS  in 
such  a  way  that  the  power  is  decreased  no  faster  than  a  certain 
power  fall  rate  {APfCSjaUrate)  to  prevent  overpressure  into  the  stack. 
Both  situations  can  be  managed  by  using  the  ESS  as  a  power  buffer. 
These  restrictions  are  formulated  in  (8).  With  regard  to  the  ESS,  the 
maximum  and  minimum  power  flows  are  also  limited,  as  is  formu¬ 
lated  in  (9).  The  maximum  power  that  the  ESS  can  supply  or  store 
depends  on  the  actual  voltage  in  the  ESS,  Vess(k),  the  maximum 
voltage  Vess.max(k),  and  the  minimum  voltage  Vess  min(k) 

Pessdwgmax  =  (VeSs(k)-V^max)Vess(k)  duringcharging  (n) 

Pess,dischmax  =  ^A) — during  discharging  (12) 
Kd 

where  Rd  is  the  ESS  internal  resistance. 

The  State  of  Energy  in  the  ESS,  SoE(k),  is  defined  as 

SoE(k):= x  ioo  (13) 

Ecap 

where  Ecap  is  the  maximum  energy  that  the  ESS  is  capable  to  store. 
Thus,  it  is  possible  to  express  the  maximum  Pess(k)  that  is  possible 
to  charge  or  discharge,  as  a  function  of  SoE(k): 


Pess,chrgmax(k )  —  kessiS0E(k)  —  SoEmax) 
Pess,dischrg  max(^)  —  kess(SoE(l<)  —  SoEmjn ) 

with 


k  Ecap 
P  d  Cr 


(14) 


(15) 


where  Rd  and  Cr  are  the  internal  resistance  and  the  capacitance  of 
the  supercapacitors,  and  k  is  a  constant  depending  in  their  partic¬ 
ular  technology.  According  to  the  sign  convention  employed,  the 
power  is  negative  when  the  ESS  is  in  charging  mode;  meanwhile  it 
is  positive  when  the  ESS  is  in  discharging  mode. 

The  internal  resistance  of  supercapacitors  is  extremely  low  and 
the  capacitance  is  exceptionally  high,  allowing  a  very  fast  opera¬ 
tion  both  during  charging  and  discharging.  In  the  supercapacitors 
employed  in  this  work,  they  are  Rd  =  0.019  £2  and  Q?  =  58F,  thus, 
kess  =  479  kW. 
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The  temporal  behaviour  of  FCS  is  fundamentally  conditioned 
by  the  compressor  dynamics  and  the  filling  dynamics.  Besides  de 
minimization  of  J,  another  important  EMS  objective  is  to  over¬ 
come  this  drawback  by  means  of  an  adequate  Energy  Management 
Strategy.  In  that  sense,  the  Energy  Management  Strategy  acts  as 
follows:  if  the  power  requested  to  the  FCS  exceeds  a  determined 
maximum  rise  rate,  then  the  power  rise  rate  that  the  FCS  actu¬ 
ally  gives  is  limited  and  the  rest  of  power  is  supplied  by  the  ESS 
if  it  is  possible  according  to  the  present  value  of  SoE(k).  If,  on  the 
contrary,  the  power  requested  to  the  FCS  exceeds  the  maximum 
fall  rate,  then  the  actual  power  fall  rate  is  limited  and  the  surplus 
power  is  absorbed  by  the  ESS  whenever  there  is  sufficient  storage 
capacity. 

3.2.  Strategy  based  on  the  FCS  efficiency  map 

One  of  the  most  relevant  characteristics  of  a  FCS  is  the  efficiency 
map  (Fig.  2),  a  chart  that  shows  how  the  efficiency  changes  with  the 
load  power.  It  is  supposed  that  the  FCS  operating  point  is  controlled 
[27]  and,  thus,  external  parameters  such  the  ambient  temperature 
have  no  influence  in  the  efficiency  map.  The  first  strategy  proposed 
in  this  work  is  a  quasi-load-following  strategy  where  the  FCS  is 
operated  in  an  advantageous  zone  where  the  efficiency  is  high. 
In  this  case,  the  operating  zone  is  constrained  between  an  infe¬ 
rior  limit  (. PfcS'io )  and  a  superior  limit  (. PfCSthi )•  The  superior  limit  is 
imposed  by  the  maximum  power  that  the  fuel  cell  can  deliver  (i.e., 
Pfcs, hi =  Pfcs, max )»  whereas  the  inferior  limit  is  determined  according 
to  the  efficiency  curve.  As  mentioned  before,  the  efficiency  of  a 
FCS  at  low  power  is  very  poor  due  to  the  parasitic  power.  Thanks 
to  the  inferior  limit,  the  Energy  Management  Strategy  avoids  this 
unfavourable  zone  (Fig.  2). 

This  EMS  gives  priority  to  the  FCS-power  supply  since  the  FCS 
is  the  primary  power  source  and  the  direct  power  flow  to  the  DC 
bus  through  the  FCS  power  converter  has  lower  losses  than  the 
indirect  way  through  the  ESS  power  converter,  the  ESS  itself  and 
again  to  the  DC  bus  through  the  power  converter  to  supply  the  load. 
Consequently,  given  the  present  PioadO<) *  if: 


PfcsJoPB  —  Pload(k)  —  Pfcs, max1! B 

(16) 

and 

A  Pfcs,  fall  rate  —  AP[oacf(k)  <  APfcsr[se  rate 

(17) 

where 

APload(k)  —  Pload(k)  ~  PloadU <  ~  1 ) 

(18) 

then,  the  FCS  operates  in  load-following  mode: 

d  (b\  PloadW 
pfcsV<)=  m 

(19) 

and 

Pess(k )  —  0 

(20) 

If,  on  the  contrary 

Pload(k)  <  Pfcs,lo{k)riB 

(21) 

then 

Pfcs(fc)  =  Pfcs,lo(k) 

(22) 

and  the  ESS  stores  the  rest  of  generated  power  if  the  ESS  is  not  too 
charged 

Pess(k)  =  —  min{|P/oad(/<)  —  Pfcsik^BldessPB/Bi  I SoE(k)  —  SoEmax\kess} 

(23) 


Fig.  3.  FCS  operating  point  for  the  strategy  based  on  efficiency  map. 


On  the  other  hand,  in  the  case  that  the  ESS  is  full  charged,  then 
PfCS(k)  =  0.  Conversely,  if 

Pload(k)  —  Pfcs,hi(k)r}B  (24) 

then 

Pfcs(k)  =  Pfcs  ,max  (25) 

and  the  ESS  supplies  the  rest  of  load  power  if  there  is  enough  energy 
in  ESS 

d  / 1  \  f  (Pload(k)  ~  Pfcs(k)riB)  rc  c(]  x  c  c  ..  \ 

PessW  —  min  <  - — — — - ,  (SoE(/<)  —  SoEm[n)keSs  f  • 

[  dB/B^ess  J 

(26) 

In  Fig.  3,  it  is  represented  how  is  determined  the  FCS  operating 
point.  Flowever,  the  transition  between  operating  points  is  limited 
by  the  maximum  rates,  thus,  finally 

p  (b)  —  f  fycs(k  —  1)  +  A Pfcs, fall  rate  AT,  if  APjcs(k)  <  APfcsjau  rate 

^  \  Pfcs ( ^  —  1 )  +  APfcs  riserate  AT,  if  APfcs(k)  >  APfCS  rise rate 

(27) 

where  APfCS  =  PfCS(1<)  -  PfCS(k  -  I),  with  PfCS(k)  as  it  was  previously 
calculated  in  (19),  (22)  or  (25),  according  to  the  conditions  (16),  (21 ) 
or  (24),  respectively.  The  power  Pess(k)  is  calculated  as  is  indicated 
in  (23)  and  (26). 

3.3.  Improved  strategy  based  on  the  FCS  efficiency  map 

The  second  strategy  based  on  the  FCS  efficiency  map  is  a  strat¬ 
egy  which  operates  the  FCS  preferably  in  its  point  of  maximum 
efficiency  in  order  to  improve  the  hydrogen  economy,  although  the 
final  operating  point  of  the  FCS  is  determined  based  on  the  actual 
power  demand  and  the  state  of  energy  of  the  ESS.  The  FCS  power 
command  is  determined  according  to  the  following  rules.  If  the  load 
power  is 

PfcsJoPB  —  Pload(k)  —  Pfcs,hi^lB  (28) 

and,  the  SoE  is 

SoEi0  <  SoE(k)  <  SoEhi  (29) 

where  Pfcs>hi  is 

Pfcs, hi  —  Pfcs,  max  PB^fcs,  hi  >  (30) 

and  Xfcs  hi  is  a  fraction  of  the  maximum  FCS  power;  then,  the  FCS  is 
operated  in  its  point  of  maximum  efficiency 

Pfcs(k)  —  Pfcs,  max  eff 


(31) 
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The  remaining  power  to  achieve  the  load  demand  flows  from 
or  to  the  ESS  according  to  (26)  if  Pi0ad(k)>Pfcs,maxeff  (discharging 
mode),  or  (23)  if  Phad(k)  <  Pfcs,max  eff  (charging  mode). 

If  the  load  power  is 


Pfcs,hi4B  —  Pload(k)  —  Pfcs^axPB 

(32) 

and  the  SoE  is 

SoEh  <  SoE(k)  <  SoEhi 

(33) 

then,  the  FCS  is  operated  in  load  following  mode 

P/cs(k)=P'°ad(k) 

(34) 

and  Pess0<)  is  as  indicated  in  (23)  or  (26). 

On  the  other  hand,  if 

Table  1 

Vehicle  model  specifications. 


Specification 

Symbol 

Value 

Vehicle  total  mass 

M 

1109  kg 

Frontal  area 

Af 

2  m2 

Drag  coefficient 

Cd 

0.335 

Coefficient  of  rolling  friction 

fr 

0.009 

Air  density 

Pa 

1.2  kg  nrr3 

Gravity 

g 

9.8  ms-2 

3.4.  Strategy  based  on  constrained  nonlinear  programming 


In  this  strategy,  an  optimization  problem  with  linear  constraints 
is  resolved  at  each  sampling  period  AT  where  a  nonlinear  cost  func¬ 
tion,  which  represents  the  hydrogen  consumption,  is  minimized. 
The  problem  can  be  put  into  the  form  [28] 


min l(x)  subject  to  H(X)  =  0  and  G(X)  <  0 
where  the  vector  X  is 


(41) 


(42) 


The  cost  function  J(X)  represents  the  hydrogen  consumption  in 
the  FCS 


J(X)  =  F(Pfcs(k))  (43) 

where 


F(Pfcs(k))  =  ConsH2(Pfcs(k))PfCS(k)AT,  (44) 

and  ConsH2(PfCS(k))  is  the  hydrogen  consumption  (gWh-1)  as  a 
function  of  Pfcs(k).  This  relationship  is  shown  in  the  form  of  a  con¬ 
sumption  map  in  Fig.  2. 

The  set  of  constraints  H(X)  =  0  represents  the  power  balance  in 
the  DC  bus  as  in  (6)  and  the  set  of  constraints  G(X)  <  0  represents  the 
limitations  in  PfCS(k)  and  Pess(k).  The  PfCS(k)  is  limited  in  its  maximum 
and  minimum  value,  and  in  the  maximum  rise  rate  and  fall  rate,  thus 


Pioad(k)  —  Pfcs,maxPB  andSoF(k)  <  SoEfc  (35) 

or 

SoE(k)  <  SoElo  (36) 

then,  the  FCS  is  operated  at  its  maximum  power 

Pfcs(fc)  —  Pfcs,  max  (37) 

and  Pess(k)  is  as  indicated  in  (26).  If,  on  the  contrary 

PioadW  <  P/cs, to Vb  and SoE(k)  >  SoElo,  (38) 

or 


Pfcs,max(k)  —  max[Pjcs  hf’,  P(k  1 )  +  APfcsjaurateAT]  (45) 

Pfcs,min(k)  =  min[P^cs  /0,  P{k  —  1)  +  APfcsriserateAT]  (46) 

On  the  other  hand,  PeSsO< )  is  limited  to  its  maximum  or  minimum 
value  depending  on  the  actual  SoE(k).  The  SoE(k)  is  limited 

SoEmin  <  SoE(k)  <  SoEmax  (47) 

thus 

kess(SoE(k)  —  SoEmax)  <  Pess(k)  —  kess{SoE(1< )  —  SoFmin)  (48) 

where  kess  is  the  constant  defined  in  (15). 


SoE(k)  >  SoEhi  (39) 

then,  the  FCS  is  operated  at  its  lower  operating  point 

Pfcs(k)  =  Pfcs,io ,  (40) 

and  Pess(k)  is  as  in  (23).  Additionally,  if  Pioad(k)  =  0  Vte  [kl,k2]  with 
(k2  -kl)>  T0ff ,  and,  SoE(k)  >  SoEhi  with  k  >  k2 ,  then,  the  FCS  is  turned 
off  to  avoid  unnecessary  hydrogen  consumption  because  the  par¬ 
asitic  losses  in  the  FCS.  Fig.  4  indicates  the  FCS  operating  point  as 
a  function  of  the  SoE(k)  and  the  load  power  Pioad0<) •  In  the  same 
way  that  in  the  previous  strategy,  the  transition  between  operating 
points  is  realized  according  to  the  constraints  concerning  the  max¬ 
imum  fall  power  rate  and  the  maximum  power  rate  as  indicated  in 
(27). 


4.  Simulation  results 

The  enunciated  strategies  were  tested  by  simulation  with  a 
model  of  a  hybrid  system  consisting  in  the  vehicle  described  in 
Table  1,  provided  with  a  FCS  and  a  supercapacitors-based  ESS  with 
the  principal  parameters  listed  in  Table  2.  The  supercapacitor  bank 
is  sized  in  order  to  accomplish  the  following  acceleration  require¬ 
ment:  the  vehicle  (with  total  mass  1380  kg  including  an  extra  cargo 
mass)  must  be  capable  to  accelerate  from  0  to  96.5  kph  in  10  s  with 
a  15-kW  FCS.  The  fuel  cell  is  sized  in  order  to  be  possible  to  fulfil  all 
the  Standard  Driving  Cycles  described  in  Section  4.1.  On  the  other 
hand,  the  strategy  parameters  P/C5./0  and  P^s, hi  are  adjusted  in  order 
to  maximize  the  fuel  economy  and  performance  according  to  the 
particular  strategy  and  cycle.  Their  values  are  shown  in  Table  3. 
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Table  2 

Parameters  of  the  powertrain  components. 


Parameter 

Symbol 

Value 

Number  of  supercapacitors 

Neap 

130 

Supercapacitor  module  mass 

Mmod 

0.408  kg 

ESS  specific  energy 

Eess,spec 

6Whkg_1 

ESS  specific  power 

P ess, spec 

2.5  kWkg-1 

Maximum  ESS  energy 

Eess.max 

318.2  Wh 

Maximum  ESS  power 

Pess,max 

132.6  kW 

Maximum  SoE 

SoEmax 

1 

Minimum  SoE 

SoEmin 

0.2 

High  limit  of  SoE 

SoEhi 

0.9 

Low  limit  of  SoE 

SoEi0 

0.3 

Maximum  FCS  power 

Pfcs.max 

15  kW 

Maximum  FCS  rise  rate  power 

APfcs,rise  rate 

600WS-1 

Maximum  FCS  fall  rate  power 

APfcsfall  rate 

-900  Ws-1 

Maximum  FCS  efficiency 

Ofcs.max 

0.6 

FCS  power  of  maximum  efficiency 

Pfcs.max  eff 

61<W 

Energy  storage  efficiency 

Oess 

0.99 

Boost  converter  efficiency 

Ob 

0.95 

Buck/Boost  converter  efficiency 

Ob/b 

0.95 

Table  3 

Parameters  PfCSj0  and  Pfcshi  used  in  the  EMSs. 

Energy  Management  Strategy 

Pfcsjc  (kW) 

PfcsM  (kW) 

Strategy  based  on  efficiency  map 

1 

15 

Improved  strategy  based  on  efficiency  map 

1 

12 

Strategy  based  on  nonlinear  programming 

2 

15 

4.1.  Testing  the  EMSs  with  standard  driving  cycles 

The  strategies  were  tested  on  four  Standard  Driving  Cycles:  the 
New  European  Driving  Cycle  (NEDC),  the  Urban  Dynamometer  Driv¬ 
ing  Schedule  (UDDS),  the  Federal  Test  Procedure  (FTP),  and  the 
Highway  Fuel  Economy  Test  (HWFET).  These  speed  profiles,  rep¬ 
resenting  urban  and  highway  scenarios,  were  originally  stated  for 
measuring  pollutant  emissions  and  gasoline  economy  of  engines 
[29]  and  are  widely  used  in  the  literature  to  evaluate  the  perfor¬ 
mance  of  EMSs  and  to  assess  the  hybridization  degree  (i.e.,  the 

parameter  that  indicates  the  relation  between  the  maximum  power 


of  the  fuel  cell  and  the  maximum  power  of  the  energy  storage  sys¬ 
tem). 

To  illustrate  the  behaviour  of  the  proposed  strategies,  the  sim¬ 
ulation  results  corresponding  to  the  NEDC  cycle  are  shown:  the 
power  split  between  the  fuel  cell  system  and  the  energy  storage  sys¬ 
tem,  and  the  evolution  of  the  SoE(t)  in  the  ESS  are  shown  in  Fig.  5  for 
the  strategy  based  on  efficiency  map,  in  Fig.  6  for  the  improved  strat¬ 
egy  based  on  efficiency  map  and  in  Fig.  7  for  the  strategy  based  on 
constrained  nonlinear  programming.  These  figures  are  only  plotted 
from  time  960  to  1180  s  for  clarity. 

4.2.  Comparative  results 

A  comparative  analysis  of  the  strategies  for  the  four  cycles,  in 
terms  of  hydrogen  consumption  per  kilometre,  is  shown  in  Fig.  8. 
The  comparative  is  done  with  respect  to  the  optimal  case  where 
the  consumption  is  minimum;  the  values  in  percentage  indicate 
the  increment  in  consumption  with  respect  to  the  optimal  case. 

The  performance  of  the  optimal  case  is  estimated  assuming  that 
the  entirely  cycle  is  known  a  priori  and  operating  the  FCS  during  the 
entirely  cycle  with  maximum  efficiency.  The  proposed  procedure 
to  evaluate  this  minimum  consumption  is  based  on  the  following 
assumptions: 

•  The  capacity  of  storage  in  the  ESS  is  sufficient  to  recover  all  the 
available  energy  from  regenerative  braking. 

•  The  friction  brake  is  not  employed  during  the  entire  cycle. 

•  The  power  consumption,  when  the  FCS  is  turned  off,  is  null. 

•  The  FCS  is  already  hot  when  the  cycle  starts. 

•  The  entire  vehicle  mass  is  fixed  for  all  the  cycles:  mveh  =  1109  kg. 

The  procedure  to  determine  the  optimal  hydrogen  consumption 
is  as  follows.  A  simulation  is  performed  for  each  driving  cycle  in 
such  a  way  that  the  FCS  works  alternately  in  two  operating  points, 
namely  “On”  and  “Off’,  according  to  the  actual  SoE(k):  (i)  when 
SoE(k)<SoE(0),  the  FCS  is  operated  at  its  point  of  maximum  effi¬ 
ciency  (the  “On”  point)  and,  (ii)  when  SoE(k)>SoE(0),  the  FCS  is 
turned  off  (the  “Off’  point).  As  a  result,  the  final  state  of  energy 
(. SoE(N ))  is  the  same  as  the  initial  state  of  energy  ( SoE(0 )).  Thus, 
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Fig.  5.  Power  split  running  on  NEDC  and  the  evolution  of  the  SoE  using  the  strategy  based  on  efficiency  map. 
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Fig.  6.  Power  split  running  on  NEDC  and  the  evolution  of  the  SoE  using  the  improved  strategy  based  on  efficiency  map. 
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Fig.  7.  Power  split  running  on  NEDC  and  the  evolution  of  the  SoE  using  the  strategy  based  on  constrained  nonlinear  programming. 


the  hydrogen  consumption  only  accounts  the  necessary  hydrogen 
to  fulfil  the  cycle  of  duration  N.  The  advantages  of  this  method  are 
considered  in  Section  6. 

In  Fig.  9,  it  is  shown  a  scheme  with  the  FCS  operation  to  perform 
the  analysis  previously  described  and  the  results  are  recollected  in 
Table  4.  Since  the  FCS  is  either  working  at  the  point  of  maximum 
efficiency  or  it  is  turned  off  (with  the  assumption  that,  when  the 
FCS  is  turned  off,  the  power  consumption  is  null),  the  minimum 
consumption  is  guaranteed.  In  addition,  in  Fig.  10  is  included  the 
performance  corresponding  to  the  pure  fuel  cell  case  where  there 
is  no  hybridization.  This  analysis  is  performed  with  a  37.5-kW  FCS, 
a  power  that  is  sufficient  to  fulfil  the  four  cycles  in  study. 


Table  4 

Driving  cycles  comparison  based  on  optimal  hydrogen  economy. 


NEDC 

UDDS 

FTP 

HWFET 

Hydrogen  consumption  (g) 

77 

80 

119 

111 

Hydrogen  consumption  (g  km-1 ) 

7.0 

6.6 

6.7 

6.7 

Braking  recovery  (Wh) 

194.2 

310.4 

452.4 

85.8 

Braking/hydrogen  ratiot  (%) 

7.6 

11.7 

11.4 

2.3 

t  The  braking/hydrogen  ratio  is  defined  as  the  ratio  between  the  energy  recovered 
from  braking  in  the  ESS  and  the  energy  of  the  Lower  Heating  Value  energy  of  the 
hydrogen  consumed  in  the  FCS. 
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■  Optimal  case  □  Pure  fuel  cell  case 

□  Strategy  based  on  efficiency  map  □  Improve  strategy  based  on  efficiency  map 

□  Strategy  based  on  constrained  nonlinear  programming _ 


NEDC  UDDS  FTP  HWFET 

Driving  cycle 


Fig.  8.  Comparative  between  the  hydrogen  consumption  of  the  proposed  strategies,  the  optimum  case  and  the  pure  fuel  cell  case.  The  percentages  in  red  indicate  the  increment 
in  consumption  with  respect  to  the  optimal  case. 


Fig.  9.  Scheme  showing  the  On/Off  operation  of  the  FCS  model  to  perform  the  analysis  of  the  minimum  hydrogen  consumption. 


4.3.  Influence  of  the  initial  SoE 

In  general,  the  final  state  of  charge  SoE(Nc)  is  different  to  the 
initial  state  of  charge  SoE(O),  resulting  in  a  gain  or  a  loss  of  energy 
in  the  ESS  over  the  driving  cycle.  Because  of  that,  the  results  are 


corrected  in  order  to  compare  the  results  correctly.  The  corrected 
consumption  of  hydrogen  is  based  on  the  assumption  that  after 
the  cycle,  the  FCS  is  continuing  running  in  its  point  of  maximum 
efficiency  until  the  SoE  reaches  again  the  initial  value.  Thus,  the 
corrected  consumption  results 


□  Strategy  based  on  efficiency  map  □  Improved  strategy  based  on  efficiency  map 

□  Strategy  based  on  constrained  nonlinear  programming 


NEDC  UDDS  FTP  HWFET 

Driving  cycle 


Fig.  10.  Hydrogen  savings  with  respect  to  the  pure  fuel  cell  case. 
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Initial  SoE  [%] 

Fig.  11.  Influence  of  the  initial  SoE  on  the  hydrogen  consumption  running  on  UDDS. 
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Fig.  12.  Influence  of  the  initial  SoE  on  the  hydrogen  consumption  running  on  HWFET. 


Fig.  13.  Scheme  showing  the  experimental  setup  to  validate  the  EMSs. 
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Cons  H2,  corrected  —  ConsH2  + 


AEess  stored 

LHVh 2  rjB  r]B /B  hess  Pfcs,  max 


(49) 


where,  Cons  H2  is  the  cumulative  consumption  of  hydrogen  over 
the  cycle  previous  to  the  correction,  AEess,  stored  is  the  difference 
of  the  energy  stored  in  the  ESS  at  the  end  of  the  cycle  with  the 
energy  at  the  beginning  of  the  cycle  (positive  if  SoE(Nc)  >  SoE(O)  and 
negative  in  the  opposite  situation),  LHVH2  is  the  low  heating  value 
of  hydrogen,  PfCS,max  is  the  maximum  efficiency  of  the  FCS  and  pess 
is  the  ESS  efficiency.  In  Figs.  11  and  12,  it  is  shown  the  influence  of 
the  initial  SoE  on  the  hydrogen  consumption  running  on  UDDS  and 
HWFET,  respectively. 


5.  Experimental  validation 

This  section  presents  the  implementation  of  the  previously 
developed  EMSs  in  an  experimental  test  setup.  The  arrangement  of 
this  Section  is  as  follows.  First,  the  experimental  setup  is  described. 
Then,  the  experimental  results  are  presented  showing  the  tempo¬ 
ral  behaviour,  the  performance  in  terms  of  hydrogen  consumption, 
and  the  comparative  between  EMSs.  Finally,  the  conclusions  of  this 
section  are  presented. 

5  A.  Description  of  the  experimental  setup 

The  experimental  setup  is  a  test  station  in  which  some  com¬ 
ponents  of  the  FCHV  are  actually  present,  some  are  emulated,  and 
some  others  are  simulated.  The  objective  of  this  setup  is  to  repro¬ 
duce,  as  close  as  possible,  the  FCHV  behaviour  in  order  to  validate 
the  EMSs.  In  this  experimental  setup  (shown  in  Fig.  13),  the  FCS  is  a 
NEXA  power  module  from  Ballard  of  1.2  kW.  The  power  module  is 
fed  with  hydrogen  from  a  pressurized  tank  and  the  hydrogen  flow 
is  measured  with  a  mass  flow  meter  (Bronkhorst,  model  F-201  AC), 
with  a  maximum  flow  capacity  of  50  SLPM. 

The  vehicle  power  consumption  is  emulated  with  a  pro¬ 
grammable  electronic  load  (Hocherl  &  Hackl,  model  ZS1606),  a 
power  device  that  can  support  up  to  1600  W  @  60  V,  150  A.  This 
load  is  commanded  by  the  Host  computer,  which  is  in  charge  of 
carrying  out  the  EMSs.  The  ESS  is  simulated  with  a  model  devel¬ 
oped  in  the  LabVIEW  environment.  LabVIEW  is  a  platform  for  a 
visual  programming  language  from  National  Instruments.  This  plat- 


Table  5 

Principal  specifications  of  BOOSTCAP  BPAK0350-15E  supercapacitors  from  Maxwell 
Technologies. 


Specification 

Value 

Capacitance,  Cr  (F) 

58 

Rated  voltage,  UR  (V) 

15 

Maximum  energy  (Wh  kg-1 ) 

3.63  (Ur  =  15  V) 

Specific  power  (W  kg-1 ) 

2850  (Ur  =  15  V) 

Internal  resistance,  Rd  (ft) 

0.019 

form  is  commonly  used  for  data  acquisition,  instrument  control, 
and  industrial  automation.  For  more  detail  see  [30]. 

The  supercapacitor  bank  is  composed  of  130  modules  and  each 
module  is  composed  of  six  serially  connected  cells.  The  main  char¬ 
acteristics  of  these  devices  are  listed  in  Table  5  and  correspond  to  a 
module  from  Maxwell  Technologies  [31  ].  In  this  experimental  stage 
the  supercapacitor  bank  is  scaled  in  the  same  proportion  as  is  the 
fuel  cell  (15-1.2  kW).  The  scaled  experimental  supercapacitor  bank 
is,  thus,  composed  of  10  modules,  the  maximum  energy  that  can  be 
stored  is  18.15  Wh,  and  the  maximum  power  that  can  be  delivered 
is  14.25  kW.  The  SoE  is  calculated  at  each  sampling  time  according 
to  this  model,  which  is  running  in  the  Host  computer. 

The  data  acquisition  and  control  system  is  composed  by  a  Host 
computer  and  a  computer  running  in  real  time,  namely  RTOS  com¬ 
puter.  The  RTOS  computer  communicates  with  the  input/output 
(I/O)  modules,  made  by  National  Instrument,  throughout  a  FPGA 
target  and  a  PCI  bus  in  a  second  computer.  The  two  computers  are 
connected  via  an  Ethernet  link.  The  Host  computer  also  allows  mon¬ 
itoring  the  evolution  of  the  system  variables  and  commanding  the 
system  through  a  graphical  interface  developed  in  LabVIEW.  The 
measured  variables  are  the  hydrogen  flow,  the  load  current  and  the 
stack  voltage.  This  variables  are  measured  every  200  ms. 

The  operation  of  the  experimental  setup  is  as  follows.  Firstly, 
the  FCS  operation  is  determined  according  to  the  EMSs  that  were 
described  in  Section  3.  To  that  end,  the  EMSs  were  programmed 
in  the  LabVIEW  environment.  These  applications  run  in  the  RTOS 
computer  and  assess  a  new  request  to  the  FCS,  PfCS(k),  every  1  s. 
Then,  the  PfCS(k)  is  required  to  the  NEXA  power  module  through 
the  electronic  load,  which  is  commanded  from  the  RTOS  computer. 
The  electronic  load  is  operated  in  power  mode,  i.e.,  the  load  sinks 
a  power  in  accordance  with  the  external-programming  signal  from 


Fig.  14.  Comparative  between  the  experimental  efficiency  map  of  NEXA  and  the  efficiency  map  of  ADVISOR  model. 
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Fig.  15.  Experimental  results  of  the  strategy  based  on  efficiency  map. 


the  RTOS  computer.  On  the  other  hand,  the  Pess(l<)  value  is  deter¬ 
mined  from  the  power  balance  in  the  DC  bus  as  in  (6)  and  is  utilized 
in  the  ESS  model,  running  in  the  Host  computer,  to  actualize  the 
actual  SoE(k). 

52.  Experimental  results 

The  previously  explained  experimental  setup  was  utilized  to  val¬ 
idate  the  enunciated  EMSs  running  on  the  same  four  driving  cycles 
in  Section  4.  In  the  simulation  stage,  it  was  utilized  a  15-kW  FCS  to 
fulfil  the  driving  cycles.  In  this  stage,  all  the  powers  are  scaled  down 
with  a  factor  of  12.5  to  fit  with  the  rated  power  of  the  NEXA  power 
module  (1.2  kW).  In  Fig.  14  is  shown  the  comparative  between 


the  experimental  efficiency  map  of  NEXA  and  the  efficiency  map 
of  the  FCS  model  in  the  ADVISOR  simulator.  The  experimental 
map  is  obtained  measuring  the  inlet  hydrogen  flow  rate  according 
to 

»•<*>- si ikr100  1501 

where  ConsH2  is  in  gkWh-1  and  LHVH2  =  120 kjg-1.  In  Figs.  15-17 
are  shown  the  experimental  results  corresponding  to  the  NEDC 
cycle,  showing  the  evolution  of  FCS  power  PfCS ,  measuring  the  stack 
voltage  and  the  stack  current;  the  evolution  of  SoE,  from  the  ESS 
model;  and  the  hydrogen  consumption,  measured  with  a  mass  flow 
meter. 
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Fig.  16.  Experimental  results  of  the  improved  strategy  based  on  efficiency  map. 
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Fig.  17.  Experimental  results  of  the  strategy  based  on  constrained  nonlinear  programming. 


In  Fig.  18,  it  is  shown  the  comparative  of  the  hydrogen  con¬ 
sumption  between  the  strategies  in  the  experimental  setup  and  the 
strategies  in  the  simulation  environment.  To  do  this  comparative 
the  experimental  results  are  multiplied  by  the  same  scale  factor 
used  to  scale  down  the  power.  From  this  comparative,  it  is  possi¬ 
ble  to  draw  some  observation.  Firstly,  the  hydrogen  consumption 
corresponding  to  the  experimental  setup  is  higher  than  the  corre¬ 
sponding  to  the  simulation  setup.  This  is  owing  to  the  fact  that  the 
efficiency  of  the  NEXA  is  lower  than  the  efficiency  of  the  FCS  model 
in  the  simulation  setup  as  can  be  seen  in  Fig.  14.  Secondly,  the  rel¬ 
ative  consumptions  in  the  experimental  setup  are  similar  to  the 
experimental  setup  with  some  differences  owing  to  discrepancies 
in  the  modelling  and  measurement  errors. 

On  the  other  hand,  in  Fig.  19  is  shown  the  hydrogen  savings 
achieved  using  the  EMSs  with  respect  to  the  pure  fuel  cell  case  with 
no  hybridization.  The  results  show  the  interest  of  hybridization, 
achieving  hydrogen  savings  up  to  33.7%. 


6.  Discussion 

The  analysis  of  the  hydrogen  economy  in  Fig.  8  shows  that  the 
three  strategies  have  a  good  performance  compared  to  the  opti¬ 
mal  case  where  the  entire  driving  cycle  is  known  a  priori,  which 
is  not  feasible  in  practice.  In  fact,  the  maximum  deviation  from  the 
optimal  case  is  9.6%,  from  the  strategy  based  on  efficiency  map  run¬ 
ning  on  NEDC  and  the  minimum  deviation  is  4.0%  in  the  improved 
strategy  based  on  efficiency  map  running  on  UDDS.  The  strategy 
based  on  constrained  nonlinear  programming  gives  the  best  per¬ 
formance  in  all  the  cases;  however,  the  performance  is  similar  to 
the  two  strategies  based  on  efficiency  map.  On  the  other  hand, 
compared  to  the  pure  fuel  cell  case  (Fig.  10),  the  results  show  con¬ 
siderable  hydrogen  savings  running  on  cycles  NEDC,  UDDS,  and 
FTP.  On  the  contrary,  running  on  FIWFET  the  savings  are  exiguous. 
This  cycle  is  a  highway  driving  cycle  and  one  of  its  characteristic  is 
that  the  average  deceleration  is  significantly  lower  (-0.22ms-2) 
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Fig.  18.  Experimental  validation  of  the  proposed  strategies. 
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□  Strategy  based  on  efficiency  map  □  Improved  strategy  based  on  efficiency  map 

□  Strategy  based  on  constrained  nonlinear  progamming 
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Fig.  19.  Hydrogen  savings  achieved  using  the  EMSs  with  respect  to  the  pure  fuel  cell  case. 


than  the  corresponding  to  the  other  cycles  (NEDC:  -0.79  ms-2; 
UDDS:  -0.58  ms-2;  FTP:  -0.58  ms-2).  These  results  suggest  that 
the  strategies  achieve  the  objectives  satisfactorily  when  a  signifi¬ 
cant  energy  from  braking  is  available. 

It  is  remarkable  that  it  is  possible  to  meet  the  load  power  in 
the  four  driving  cycles  with  a  15-kW  fuel  cell,  what  is  significantly 
lower  to  the  37.5-kW  fuel  cell  that  is  necessary  if  no  hybridization 
is  present.  This  advantage  would  be  translated  in  a  reduction  in  the 
production  costs  and  hydrogen  consumption. 

The  analysis  of  the  influence  of  the  initial  SoE  on  the  hydrogen 
consumption  (Figs.  11  and  12)  shows  that  the  first  strategy  and  the 
third  strategy  have  a  minimal  dependence  on  this  parameter,  which 
denotes  robustness  of  the  strategies.  On  the  other  hand,  the  per¬ 
formance  of  the  second  strategy  is  dependent  on  the  initial  SoE, 
especially  in  the  HWFET  cycle.  In  this  case,  the  performance  is  worse 
than  in  the  first  strategy  if  the  initial  SoE  is  lower  than  60%. 

The  supercapacitor  bank  is  designed  in  order  to  fulfill  an  acceler¬ 
ation  requirement  as  mentioned  in  Section  4.  However,  it  is  useful 
to  analyze  how  the  strategies  work  with  different  sizes  of  superca¬ 
pacitor  bank.  When  the  number  of  supercapacitors  is  increased,  the 
hydrogen  consumption  slightly  increases  because  of  the  increment 
in  the  vehicle  mass  is  the  preponderant  factor.  On  the  contrary, 
when  the  number  of  supercapacitors  starts  to  be  decreased,  the 
hydrogen  consumption  also  decreases  because  of  the  mass  reduc¬ 
tion  although  the  acceleration  test  is  not  always  fulfilled. 

Concerning  to  the  analysis  of  the  optimal  hydrogen  consump¬ 
tion,  in  some  works,  e.g.,  [19,32-33],  it  is  performed  through  the 
use  of  the  Dynamic  Programming  (DP)  technique,  a  procedure  that 
allows  to  determine  a  global  optimal  operation  of  the  system  for 
a  given  load  profile  by  evaluating  all  possible  control  sequences 
in  a  systematic  way.  However,  a  disadvantage  of  this  technique  is 
the  relative  long  computation  time  due  to  the  large  required  grid 
density.  The  grid  is  the  result  of  the  discretization  of  time  and  the 
state  variables  in  appropriate  grid-levels.  Thus,  in  order  to  obtain 
accurate  results  a  high  density  grid  has  to  be  considered.  On  the 
contrary,  in  our  work,  this  issue  is  resolved  without  discretization 
in  the  state  variables,  reducing  the  computation  time  without  losing 
accuracy. 

7.  Conclusions 

In  this  work,  Energy  Management  Strategies  for  Fuel  Cell  Hybrid 
Vehicles  were  addressed.  The  objectives  of  the  EMSs  were  enun¬ 
ciated  and  three  strategies  were  proposed.  These  strategies  are 


based  on  the  knowledge  of  the  efficiency  map.  The  first  strategy 
is  a  quasi-load-following  strategy  in  which  the  FCS  is  operated  in 
an  advantageous  zone  where  the  efficiency  is  high  whereas  the 
second  strategy  operates  the  FCS  preferably  in  its  point  of  maxi¬ 
mum  efficiency  in  order  to  improve  the  hydrogen  economy.  In  the 
third  strategy,  an  optimization  problem  with  linear  constraints  is 
resolved  at  each  sampling  period  where  a  nonlinear  cost  function, 
which  represents  the  hydrogen  consumption,  is  minimized. 

Firstly,  the  EMSs  were  evaluated  in  a  simulation  environment 
using  four  standard  driving  cycles  and,  then,  the  EMSs  were  tested 
in  an  experimental  setup  reproducing  the  FCVH  behaviour  were 
some  components  are  actually  present,  some  are  emulated  and 
some  others  are  simulated. 

The  results  show  that  using  the  proposed  EMSs  it  is  possible  to 
achieve  a  high  reduction  in  the  hydrogen  consumption  in  the  range 
of  17%  average.  This  is  achieved  operating  the  FCS  properly  and  tak¬ 
ing  advantage  of  the  energy  from  regenerative  braking.  To  evaluate 
the  savings,  the  hydrogen  consumptions  are  compared  with  two 
references:  the  pure  fuel  cell  case  with  no  hybridization  and  the 
optimal  case  with  minimum  consumption  where  the  driving  cycle 
is  known  a  priori. 
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